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Abstract 

We present preliminary results of a search for charmless two-body B decays to charged pions and 
kaons using data collected by the BABAR detector at the Stanford Linear Accelerator Center's PEP- 
II storage ring. In a sample of 8.8 million produced BB pairs we measure the branching fractions 
B{B° -> vr+Tr-) = (9.3±|6+i-2) x 10 -6 and #(#0 ^ k+tt') = (12.5^;°!^) x lO" 6 , where the 
first uncertainty is statistical and the second is systematic. For the decay B° — > K + K~ we find no 
significant signal and set an upper limit of B(B° — ► K + K~) < 6.6 x 10~ 6 at the 90% confidence 
level. 
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1 Introduction 



Measurements of the branching fractions for the rare charmless decays B° — ► h + h~ (h = ir,K)^ 
provide important information in the study of charge-parity (CP) violation. In principle, the 7r + 7r~ 
decay mode can be used to extract the angle a of the Unitarity Triangle through the phenomenon 
of B°-B° mixing. However, in addition to the dominant b — > uW~ tree amplitude, this decay 
includes the b — > dg penguin amplitude and the determination of a is subject to large theoretical 
uncertainties if the penguin contribution is non-negligible pj. In the presence of significant "penguin 
pollution," additional measurements of the isospin-related decays B~ — ► 7r~7r and B — > 
provide a means of measuring a cleanly in the tttt channel H . 

The decay B°^K+ 7r is dominated by the b — > sg penguin amplitude and provides an estimate 
of the scale of penguin pollution in the 7r + 7r _ decay. Recent results from the CLEO collaboration 
indicate that the decay rate for B°^K + tt~ is significantly larger than the rate for B°— >ir + ir~ ||, 
implying that the penguin contribution in the tttt decay is indeed significant and an isospin analysis 
will be necessary to measure a accurately. The apparent enhancement of the penguin amplitude 
improves the prospects for observing direct CP violation as an asymmetry in the decay rates 
for B°^>K + 7r~ and B°^>K~it + . Precise measurement of the decay rates for tttt and Ktt decays 
is therefore of central importance. This paper describes preliminary measurements of branching 
fractions for the decays B ^7r + n~ , K + ir~ , and K + K~ with the first data collected by the BABAR 
experiment. 

2 Data sample, BABAR detector, and event selection 

The dataset used in this analysis consists of 8.9 fb -1 collected with the BABAR detector at the 
Stanford Linear Accelerator Center's PEP-II storage ring between January and June 2000. The 
PEP-II facility operates nominally at the T(4S) resonance, providing asymmetric collisions of 
9.0 GeV electrons on 3.1 GeV positrons. The dataset includes 7.7 fb -1 collected in this configuration 
(on-resonance) and 1.2 fb -1 collected below the BB threshold (off-resonance) that are used for 
continuum background studies. The on-resonance sample corresponds to 8.8 million produced BB 
events. 

The asymmetric beam configuration in the laboratory frame provides a boost (/3j = 0.56) to the 
T(45), allowing separation of the B and B decay products for time-dependent CP- violation studies. 
For the analysis described in this paper, the significant effect of the boost relative to symmetric 
collider experiments is to increase the momentum range of the B decay products from a narrow 
distribution centered at ~ 2.6GeV/c, to a broad distribution extending from 1.5 to 4.5GeV/c. 
Wherever necessary, kinematic quantities evaluated in the T(4J5) center-of-mass (CM) frame are 
denoted with an additional superscript asterisk in order to distinguish them from the corresponding 
quantities evaluated in the laboratory frame. 

BABAR is a solenoidal detector optimized for the asymmetric beam configuration at PEP-II and 
is described in detail elsewhere Q]. Charged particle (track) momenta are measured in a tracking 
system consisting of a 5-layer, double-sided, silicon vertex tracker and a 40-layer drift chamber 
filled with a gas mixture of helium and isobutane, both operating within a 1.5 T superconducting 
solenoidal magnet. Photons are detected in an electromagnetic calorimeter consisting of 6580 
CsI(Tl) crystals arranged in barrel and forward endcap subdetectors that also operates within the 
magnetic field. The iron used for the magnet yoke is segmented and instrumented with resistive 

1 Charge conjugate decay modes are assumed throughout this paper. 
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plate chambers, providing muon identification and (in conjunction with the calorimeter) neutral 
hadron detection. 

In this analysis, tracks are identified as pions or kaons by the Cerenkov angle 9 C measured by 
a Detector of Internally Reflected Cerenkov light (DIRC). The DIRC system is a unique type of 
Cerenkov detector that relies on total internal reflection within the radiator to deliver the Cerenkov 
light outside the tracking and magnetic volumes. The radiator consists of 144 axially aligned 
synthetic quartz bars located just inside the inner radius of the calorimeter. The bars extend 
outside the solenoid flux return in the backward direction, where the Cerenkov ring is imaged by 
an array of ~ 11000 photomultiplier tubes. 

Hadronic events are selected based on track multiplicity and event topology. Tracks with trans- 
verse momentum greater than 100 MeV/ c are required to pass efficient quality cuts, including num- 
ber of drift chamber hit layers used in the track fit and impact parameter in the r— <j> and r-z planes, 
where the cylindrical coordinate z is aligned along the detector axis in the electron beam direction. 
At least three tracks must pass the above selection. To reduce contamination from Bhabha and 
events the ratio of second to zeroth Fox- Wolfram moments ||, R2 = H2/H0, is required 
to be less than 0.95. Residual background from tau hadronic decays is reduced by requiring the 
sphericity || of the event to be greater than 0.01. The efficiency of the event selection is domi- 
nated by the acceptance and efficiency of the track requirement, and is determined to be 70% from 
a detailed Monte Carlo simulation based on Geant321 0]. 

3 Selection of B — > ft + fa - candidates 

One of the advantages of measuring B decay parameters at the T(4S) resonance is the kinematic 
constraint provided by the initial state, where energy conservation determines that the energies of 
the B mesons in the CM frame are equal to \fs/2, where -^/s is the total e + e~ CM energy. We 
exploit this constraint by calculating an energy-substituted mass tties, where \fs/2 is substituted 
for the B candidate energy, and by calculating the energy difference AE between the B candidate 
and \/s/2 in the CM frame. 

We define tties = J (\/~s/2) 2 — p* B , where p* B is the B candidate momentum evaluated in the 
CM frame. Because p* B is relatively small (~ 300MeV/c), the resolution on toes is dominated 
by the uncertainty in -y/i, which in turn is determined by the beam energy spread and width of 
the T(4S) resonance. Substitution of the beam energy reduces the mass resolution by one order 
of magnitude compared to the invariant mass. The mean value of toes an d its Gaussian width 
tr(mEs) are determined from a large sample of fully reconstructed B decays. We find toes = 
(5.2800 ± 0.0005) GeV/c 2 and a(m E s) = (2.6 ± 0.1) MeV/ c 2 , respectively. Our initial selection 
requires 5.22 < toes < 5.30GeV/c 2 . 

We define AE = E B — y^/2, where E* B is the B candidate energy in the CM frame. Signal 
events are Gaussianly distributed in AE with a mean near zero, while the continuum background 
falls linearly over the region of interest. For this analysis, the pion mass is assigned to all tracks and 
the Kit and KK decays have AE shifted from zero by an amount depending on the momentum 
of the tracks. From Monte Carlo simulation we find shifts of —45 and —91 MeV for the Kir and 
KK decays, respectively. The resolution on AE is estimated to be 27 ± 5 MeV based on Monte 
Carlo simulated tttt decays and the observed difference in widths between data and Monte Carlo 
B-^D°ir- decays. We require \ AE\ < 0.420 GeV. 
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4 Background suppression 



Due to the relatively small CM momenta of decay products produced from the quark transition 
b — > c, B decays to final states involving charm mesons are not a significant background to charmless 
two-body decays. After hadronic selection, the background is dominated by continuum production 
of light quarks, e + e~ — > qq(q = u,d,s,c). In the CM frame, the continuum background typically 
exhibits a two-jet structure that can produce two high momentum back-to-back tracks with an 
invariant mass near the B mass. In contrast, the low momentum of B mesons in the decay 
Y{AS)^>BB leads to a more spherically symmetric event. This topology difference is exploited 
by constructing the angle 9s between the sphericity axes, evaluated in the CM frame, of the B 
candidate and the remaining charged and neutral particles in the event. The absolute value of the 
cosine of this angle is strongly peaked near 1 for continuum events and is approximately uniform 
for BB events. We require |cos%| < 0.9, which is 87% efficient for signal events and rejects 66% 
of the continuum background. 

Further separation power between signal and continuum background is provided by a Fisher 
discriminant technique |||. The Fisher discriminant T is calculated from a linear combination of 
discriminating variables X{, 

9 

J 7 = ^2a i x i , (1) 

i=l 

where the coefficients a» are chosen to maximize the statistical separation between signal and 
background events. The nine discriminating variables are constructed from the scalar sum of the 
momenta of all charged and neutral particles (excluding the candidate daughter tracks) flowing 
into nine concentric cones centered on the i3-candidate thrust axis in the CM frame. Each cone 
subtends an angle of 10° and is folded to combine the forward and backward intervals. More energy 
will be found in the cones nearer the candidate thrust axis in jet-like continuum background events 
than in the more isotropic BB events. 

Large samples of signal and background Monte Carlo simulated events reconstructed in the titt 
mode are used to determine the Fisher coefficients. Figure |l| shows the resulting T distributions 
for signal nir Monte Carlo compared to a sample of B~^D°tt~ decays reconstructed in data, and 
continuum background Monte Carlo compared to off-resonance data. The J- distributions for both 
signal and background are parameterized by the sum of two Gaussians with separate means and 
widths. The Gaussian fits are performed on Monte Carlo samples that are independent of the 
samples used to determine the Fisher coefficients, and the same parameterization is used for all 
three signal modes. 

5 Particle identification 

Two complementary methods of exploiting the particle identification capabilities of the DIRC are 
described in this paper. The first method uses measurements of C to derive particle selectors that 
are used to identify pions and kaons on a per-track basis. The second method uses likelihood 
functions derived from 9 C measurements directly in a maximum-likelihood fit to extract the relative 
amount of each decay mode on a statistical basis. After the selection criteria described above, the 
combined acceptance and efficiency of requiring a 6 C measurement for both tracks is 76%. 

A control sample consisting of 18141 ± 140 D°^K~tt + candidate decays is used to parame- 
terize and assess the performance of both particle identification methods. A 96% pure D° sample 
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-1 1 

Fisher Discriminant 



Figure 1: The T distribution for tttt signal Monte Carlo events (solid histogram and fitted curve) 
compared to data B~ — > D°tt~ decays (filled squares), and continuum background Monte Carlo 
(dashed histogram and fitted curve) compared to off-resonance data (open squares). The Monte 
Carlo samples are independent of the samples used to train the Fisher discriminant. 



is obtained through the decay D* + — > D°irf, where the slow pion 7r+ tags the charge of the pion 
from the D° decay. Requiring a tight window around the D* + -D° mass difference minimizes con- 
tamination from incorrectly reconstructed D° mesons, where the two daughter tracks are assigned 
the incorrect particle hypotheses. The momentum of kaons and pions in the control sample ranges 
from 1.75-4.0 GeV/c, which covers 90% of the momentum range for charmless two-body B decays. 

The performance of the DIRC is summarized in Fig. ^, where we show plots of (a) C vs. p 
for tracks in the control sample and (b) the statistical separation between pions and kaons as a 
function of momentum. The separation, defined as [{6 c (n)) — (6 C (K))] /(crg c ), where {0 c {h)) are the 
Gaussian means for pion and kaon tracks and (ctq c ) is the average width, varies from 8 at 2GeV/c 
to 2.5 at 4GeV/c. Note that the quoted 6 C separation is for a single track and an average over the 
polar angle of each track is implicit in Fig. ||(b). 

The selector method of particle identification attempts to identify pions and kaons on a per-track 
basis by cutting on likelihood functions. We construct the likelihood for a given mass hypothesis 
from the Gaussianly-distributed C measurements and the Poissonian probability for the number of 
measured Cerenkov photons compared to its expectation value. Tracks with momentum below the 
kaon threshold or within 2a(9 c ) of the expected value for a proton are explicitly removed, where 
the measurement error cr(6 c ) varies with momentum and the number of photons used in the fit. 
We calculate the efficiency for one mode to be identified as another using the efficiency and mis- 
identification probabilities for pions and kaons measured in the D° control sample. The resulting 
efficiency matrix is shown in Table |]. The matrix includes the probability that a pion or kaon is 
mis-measured and removed by the proton rejection requirement. 

The global likelihood method incorporates the particle identification probabilities for pions 
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Figure 2: (a) The Cerenkov angle and (b) K—ir separation as functions of momentum for single 
tracks in the D° control sample. The separation is an average over all polar angles. 



Table 1: Particle selector efficiency and mis-identification probability for each decay mode. Errors 
are statistical only. 

Probability to be identified as 
Mode 7T7T Ktt KK 

vrvr 0.853 ± 0.005 0.109 ±0.008 0.0029 ± 0.0002 
Ktt 0.121 ±0.004 0.775 ± 0.006 0.051 ± 0.004 
KK 0.0112 ±0.0006 0.231 ± 0.008 0.704 ± 0.007 



and kaons directly in an unbinned maximum likelihood fit. The probability density functions are 
constructed from 6 C information alone, where Gaussian fits are performed to the distribution of 
measured — expected Cerenkov angle in bins of momentum in the D° control sample. To improve 
the resolution and minimize the non-Gaussian tails of the 9 C distribution, we require a minimum 
number of observed Cerenkov photons above the expected background. From the control sample we 
determine the efficiency of this requirement to be 77%, taking into account the angular correlation 
between the two daughter tracks from a B decay. The proton rejection cut applied in the selector 
method is also required for the likelihood analysis. After the minimum photon cut, the proton 
requirement is 98% efficient for h + h~ events. 

6 Analysis 

In this section we present the results of a simple cut-based analysis to determine the h + h~ yield, 
followed by two complementary methods for determining the irir, Kir, and KK yields in our data. 




2 3 4 
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In the first method we apply background suppression and particle identification cuts to isolate 
samples of events that are consistent with the tttt, Kit, or KK hypotheses. Signal yields are then 
obtained from an unbinned maximum likelihood fit to tties- I n the second method we perform 
a global maximum likelihood fit incorporating toeS) AE, and J 7 , as well as the 9 C probability 
density functions described in the previous section, to determine signal yields in all three modes 
simultaneously. The cut-based method is more transparent while the global fit has higher efficiency 
and statistical significance. In Sec. |7] we calculate branching fractions using the global fit results. 



6.1 Cut-based analysis 

In addition to the selection criteria described in Sees. || and |4], we tighten the 9 s requirement, 
|cos#s| < 0.7, and also require |cos#s| < 0.8, where 9b is the CM angle of the B candidate with 
respect to the beam. The variable cos 9b is uniform for the continuum background and follows 
a 1 — cos 2 9b distribution for signal B decays. We determine an optimal selection requirement 
of T > 0.37 by maximizing the statistical significance of the expected signal yield in a sample of 
Monte Carlo simulated signal and background events. The relative efficiency of these additional 
cuts is 44%. Within the AE-m^s plane we define a signal region \AE\ < 0.140 GeV and mEs = 
5.2800 ± 0.0052 GeV/ c 2 , and sideband regions \AE\ > 0.140 GeV and m ES < 5.27 GeV/ c 2 . The 
signal region in AE is designed to minimize contamination from three-body charmless B decays 
where one of the decay products has low momentum. 

Before applying the particle selector we first demonstrate the presence of h + h~ decays in the 
signal region. The signal yield is obtained from an unbinned maximum likelihood fit to the m-ES dis- 
tribution. The fit includes candidates passing all cuts except the requirement that the tracks have 
an associated 9 C measurement. The background shape in tties is parameterized by the empirical 
formula ||], 

f(m E s) oc m E svl - x 2 exp -£(1 - x 2 ) , (2) 

where x = 2mEs/\/s and the parameter £ is determined from a fit. Figure || shows the results of 
fitting Eq. || to the on-resonance AE sideband region, where we find £ = 22.0 ±0.5. We then fit the 
mES distribution in the AE signal region to a Gaussian with mean and width fixed to 5.280 GeV/c 2 
and 2.6MeV/c 2 , respectively, and the background shape fixed to £ = 22; only the normalizations 
are free parameters. The result is shown in Fig. |j. The fitted number of h + h~ candidates is 67 ± 11, 
where the error is the statistical uncertainty from the fit. Correcting for the total efficiency of the 
selection criteria (27%) and normalizing to the total number of BB pairs, we determine a branching 
fraction B(B° — > h + h~) = (28 ± 5) x 10~ 6 , where the error is statistical only and we assume equal 
Y(AS) branching fractions to charged and neutral B mesons. 

The presence of kaons in the signal region can be demonstrated by plotting the variable 
(6> c — c (ir)) j (9 c (tt) — 9 C (K)), where 9 c (ir) and 9 C (K) are the expected values. This variable peaks 
at (— 1) for real pion (kaon) tracks. In Fig. |5] we show the distribution of this variable for B- 
candidate tracks in the mES signal region before and after subtracting the distribution obtained 
from the toes sideband region. A clear kaon peak remains after subtraction. 

We now decompose the signal into its constituent components by using the particle selector to 
separate the sample into three subsamples corresponding to the different channels. We perform 
three separate fits similar to the h + h~ fit. After correcting for the efficiency and mis-identification 
probabilities in Table [l], we find 25 ± 8 tttt, 26 ± 8 Kir, and 1.2l?j KK decays. In Fig. § we show 
the fit results for mES an d the AE distribution in the mES signal region for all three modes. 
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Figure 3: Result of fitting of Eq. || to the AE sideband region in on-resonance data. The result is 
f = 22.0 ±0.5. 
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Figure 4: Result of fitting the wes signal region in the h + h~ sample to the sum of Eq. |], with 
£ = 22, and a Gaussian signal peak with toes = 5.28 GeV/ c 2 and cr(mEs) = 2.6 MeV/ c 2 . The dashed 
curve shows the background parameterization only. 



To facilitate comparison with the global maximum likelihood fit, we determine the non-common 
systematic errors for the cut-based analysis. We vary the mean and width of the mES signal 
and background distributions within their estimated uncertainties. Systematic uncertainty on the 
particle identification method is estimated by comparing with various selector definitions, including 
the probability density functions used in the global fit (below). Uncertainty in the T shape is 
estimated by varying the cut and by substituting the shape obtained from the B — > D ir~ data 
sample (Fig. [l]). Table || summarizes these uncertainties. 
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Figure 5: The measured Cerenkov angle minus its expectation value for the pion hypothesis, 
normalized to the expected separation between pions and kaons. True pions (kaons) appear as 
a peak at (—1). (a) shows the distribution for ^candidate tracks in the signal region of AE and 
the signal (solid) and sideband (dashed) regions of toes; while (b) shows the signal region after 
sideband subtraction. Note that both tracks from each B candidate are included in these plots. 
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Figure 6: Fit results for ?h,es in the (a) tttt, (b) Kn, and (c) KK samples defined by the particle 
identification selector described in Sec. |5|. Also shown are the AE distributions for the (d) tttt, (e) 
Ktt, and (f) KK samples in the mss signal region. 
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Table 2: Systematic errors (%) in the cut-based analysis that are not common to the global fit 
analysis. "PID" refers to the method of particle identification. 



Mode 


Bkg m ES 




0"O™Es) 




T cut 


PID 


Total 




0.4 


+U.4 
-1.2 


0.8 


1.9 


11 


9 


±14 


K+ir- 


0.4 


+1.7 
-3.0 


0.8 


1.9 


8 


9 


±13 


K+K- 
h + h- 


0.5 
0.5 


+23 
-27 
+1.2 
-2.7 


+2.7 
-5.3 
+1.3 
-1.2 


1.9 
1.9 


14 

6 


14 


+31 
-34 

7 



6.2 Global fit 

We perform an unbinned maximum likelihood fit using m-Es, AE, T , 9±(pi), and 02(P2), where 0\ 
and 02 are the Cerenkov angles for each track and p\ and p2 are the momenta. The likelihood C is 
defined as 

N 

C = e- N ']]Vi(m ES ,AE,f,6 1 (p 1 ),62(p2)\N n ^N Kn ,N nK ,N KK ,N hkg ), (3) 
i=i 

where N^ kg is the number of continuum background events and V% is the probability for the ith 
candidate assuming the total yield 

N' = N W7r + N Kjr + N nK + N KK + N hkg . (4) 

In the Ktt mode we fit separately for the two possible combinations (irK or Kit). The term (e~ N ') 
derives from the Poissonian probability of observing N total events when N' are expected. The 
probability for a given candidate is the sum of the signal and background terms 

V i {m ES ,AE,F,e 1 {p 1 ),62(p 2 )\N 7TW ,N Kn ,N nK ,N KK ,N hkg ) =£^7^, (5) 

k 

where the index k represents the five fit components and V\ is the product of probability density 
functions for uies, AE, J 7 , and 62^2) ■ 

The fit includes all candidates satisfying the selection criteria described in Sees. Q |3|, and 
as well as the requirement on the number of Cerenkov photons above background and the proton 
rejection cut described in Sec. ||. Due to the use of the discriminating variable AE, the global fit 
is much less susceptible to contamination from three-body B decays. The signal region is therefore 
expanded to —0.200 < AE < 0.140 GeV. The quantity — log£ is minimized with respect to the fit 
parameters. The resulting signal yields are 29l^ irir, 38^g Ktt, and 7^ KK decays. As a visual 
cross check, in Fig. [7| we plot the ttt-esj AE, and C distributions after the additional requirements 
on cos 6s, cosOb, and T applied in the cut-based analysis, and overlay the global fit results after 
rescaling by the relative efficiency for these cuts. 

Systematic errors on the fit results are estimated by varying the signal and background prob- 
ability density functions for ttiesi AE, and 8 C within their errors. The AE width is significantly 
different between data and Monte Carlo simulated B~ — > D°tt~ decays and there is evidence for a 
few MeV shift of the mean value in the negative direction. To be conservative we vary the width by 
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Figure 7: The distributions of (a) toes, (b) AE, and (c) C (solid histograms) after cutting on 
cos 6 Si cos 9b, and T. The global fit results are overlaid after scaling by the relative efficiency of 
the additional cuts. The shaded region is the total signal contribution, while the dashed and dotted 
lines are the nir and Kir components, respectively. 



Table 3: Systematic errors (%) on the signal yield from the global fit that are not common to the 
cut-based analysis. Uncertainty on 9 C includes mean and width. 



Mode 


Bkg to es 


Bkg AE 


(m E s) 


0"(mEs) 


(AE) 


a(AE) 


^0D U 7T) 


e c 


Total 


7T + 7T~ 


0.2 


0.4 


+U.Y 

-2.0 


0.3 


+1.9 


+7 

-10 


7.0 


+3.5 
-2.4 


+11 

-13 




0.2 


0.2 


+1.0 

-5 


+0.3 
-0.4 


-0.7 


+5 

-9 


5 


+1.0 
-1.4 


+7 
-12 


K+K- 


0.3 


6 


+1.3 
-3.1 


+2.1 
-2.5 


-3.1 


+2.7 
-2.9 


22 


2.0 


+23 
-24 
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±5MeV and refit assuming {AE) = — 5MeV. The uncertainty due to the shape of T is determined 
by using the shape obtained from data B — * D°tt~ decays. Table [3| summarizes the systematic 
errors that are not common to the cut-based analysis. 

6.3 Comparison and cross-checks 

In order to compare results for the two methods, we correct the signal yields obtained in the cut- 
based analysis by the relative efficiency of the |cos#s| < 0.7, |cos6*b| < 0.8, and T > 0.37 cuts 
(44%), and the global likelihood fit yields by the efficiency of the Cerenkov photon cut and proton 
rejection cuts (76%). In order to compare total yields to the h + h~ fit result we also correct by the 
common requirement that both tracks have an associated 9 C measurement (76%). The resulting 
comparison is summarized in Table ^[ Given that the two samples are not 100% statistically 
correlated, we find the agreement to be satisfactory. 

As a cross-check we perform the global fit after applying the particle selector and the additional 
selection criteria used in the cut-based analysis (including |Ai5| < 0.140 GeV). The AE and 9 C 
distributions are used in the fit, the T distribution is not. The results, summarized in Table || 
indicate that there is no significant migration between signal categories. This result confirms the 
consistency in the use of 6 C in both analyses. 

Table 4: Comparison of signal yields corrected for the relative efficiency of cut-based and global 
fit analyses. The first uncertainty is statistical, the second is systematic. For total h + h~ yield we 
compare the sum of central values for the individual modes in the two analyses with the h + h~ fit 
result (Fig. |). 



Mode 


cut-based 


global fit 


h + h~ fit 


7T + 7T~ 


76 ± 23 ± 11 


ou -12-7 




K+ir- 


77 ± 25 ± 10 


R7 +16+5 
0'_l4-8 




K + K- 


4^±1 


13±|±3 




h+h- 


156 


129 


151 ±24 ±11 



Table 5: Signal yields from global fits to the particle-selected samples. Note that these results have 
not been corrected for the 76% efficiency of the cut on minimum number of Cerenkov photons. 







Selected as 




Mode 


TTTT 


Kit 


KK 


TTTT 


21 ±5 


u - u -o.o 


o +u - b 

u - u -o.o 


Ktt 


o o+ - 7 
u - u -o.o 


24 ±5 


o o +0 - 7 
u - u -o.o 


KK 


n n +0 - 5 
u - u -o.o 


1.7 ±2.3 


2.1 ± 1.9 
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Table 6: Summary of branching fraction results for the global likelihood fit. Shown are the central 
fit values Ns, the statistical significance, and the measured branching fractions B. For the KK 
mode, the 90% confidence level upper limit on the signal yield is given in parenthesis. There is a 
common efficiency of 0.35 db 0.02 for all three modes. For N$ and B the first error is statistical and 
the second is systematic. Charge conjugate modes are assumed. 



Mode 


N s 


Stat. Sig. 


(a) B(10-») 


7T + 7T~ 




5.7 


y -°-2.3-1.4 


K+TT- 


38±1±1 


6.7 


in tc+3.0+1.3 
iz -' J -2.6-1.7 


K + K- 


7±l (< 15) 


2.1 


< 6.6 



7 Determination of branching fractions 

We determine branching fractions for ir + ir~ and K + tt~ decays and an upper limit for the K + K~ 
decay using the results of the global likelihood fit. The individual efficiencies were reported in 
previous sections. The total efficiency is 0.35 ± 0.02, where the error is combined statistical and 
systematic. Branching fractions are calculated as 

3 = 7ih> (6) 

6 BB 

where Ns is the central value from the fit, e is the total efficiency, and N B g is the total number 
of BB pairs in our dataset. Implicit in Eq. |6] is the assumption of equal branching fractions for 
T(AS) — > B°B° and T(4S) — * B + B~. The results are summarized in Table ||. In addition to 
the systematic uncertainties listed in Table [31 the total error includes uncertainty on the tracking 
efficiency (2.5% per track) Q, the shape of cos 9s (3%), and the number of BB events (3.6%). 

The statistical significance of a given signal yield is determined by setting the yield to zero and 
maximizing the likelihood with respect to the remaining variables. The results are give in Table §. 
Fig. [H shows the na likelihood contour curves, where a represents the statistical uncertainty only. 
The curves are computed by maximizing the likelihood with respect to the remaining variables in 
the fit. For the KK mode we calculate the 90% confidence level upper limit yield and decrease the 
efficiency by the total systematic error (24%) before calculating the upper limit branching fraction. 



8 Summary 

We have performed a search for charmless two-body B decays to charged pions and kaons. The 
statistical significance of the signal yields are 5.7, 6.7, and 2.1 standard deviations for the tttt, Kit, 
and KK decay modes, respectively. For the decay modes B° — ► ir + ir~ and B° — ► K + ir~ we measure 
preliminary branching fractions of (O.Sj^'sii 4) x 10 -6 and (12.5jT 2 6— 1 7) x 10 -6 > respectively, where 
the first uncertainty is statistical and the second is systematic. Since the KK yield is not significant 
we calculate the 90% confidence limit and find B(B° — > K + K~) < 6.6 x 10 -6 . 
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Figure 8: The central value (filled circle) for B(B — ► ir + ir~ -) and B(B° -> K+vr") along with the 
nff statistical contour curves for the global likelihood fit, where no corresponds to a change of n 2 
in — 21og£. 
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